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Abstract 

As discussed recently by Hooper and Tait, the singlino-like dark matter in the Minimal Su- 
persymmetric Standard Model (MSSM) extended by a singlet Higgs superfield can give a perfect 
explanation for both the relic density and the Pamela result through the Sommerfeld-enhanced 
annihilation into singlet Higgs bosons (a or h followed by h — aa) with a being light enough to 
decay dominantly to muons or electrons. In this work we analyze the parameter space required 
by such a dark matter explanation and also consider the constraints from the LEP experiments. 
We find that although the light singlet Higgs bosons have small mixings with the Higgs doublets 
in the allowed parameter space, their couplings with the SM-like Higgs boson hsm (the lightest 
doublet-dominant Higgs boson) can be enhanced by the soft parameter A,, and, in order to meet 
the stringent LEP constraints, the hsm tends to decay into the singlet Higgs pairs aa or hh in- 
stead of bb. So the hsm produced at the LHC will give a multi-muon signal, hgm —> aa > 4p or 


hsm > hh —> 4a — 8u. 


PACS numbers: 14.80.Ly,11.30.Pb,95.35.+d 


I. INTRODUCTION 


The experiment Pamela has observed an excess of the cosmic ray positron in the energy 
range 10-100 GeV [1], which is hard to be explained by the conventional cosmic ray source 
[2]. While there may exist some mundane explanations like pulsars [3] and the acceleration 
of positron secondaries in cosmic ray acceleration regions [4], the dark matter interpretation 
[5, 6] is especially interesting since it may be related to new physics to be probed at the 
LHC. 

To explain the Pamela excess by the dark matter annihilations, there are some challenges. 
First, the dark matter must annihilate dominantly into leptons since Pamela has observed 
no excess of anti-protons [1] (However, as pointed in |7], this statement may be not so solid 
due to the significant astrophysical uncertainties associated with their propagation). Second, 
the explanation of Pamela excess requires an annihilation rate which is too large to explain 
the relic abundance if the dark matter is produced thermally in the early universe. To tackle 
these difficulties, a new theory of dark matter was proposed in [6]. In this new theory the 
Sommerfeld effect of a new force in the dark sector can greatly enhance the annihilation 
rate when the velocity of dark matter is much smaller than the velocity at freeze-out in the 
early universe, and the dark matter annihilates into light particles which are kinematically 
allowed to decay to muons or electrons. 

The above fancy idea is hard to realize in the popular Minimal Supersymmetric Standard 
Model (MSSM) because there is not a new force in the neutralino dark matter sector to 
induce the Sommerfeld enhancement and the neutralino drak matter annihilates largely to 
final states consisting of heavy quarks or gauge and/or Higgs bosons [8, 9]. However, as 
discussed in [10], in the extension of the MSSM by introducing a singlet Higgs superfield, 
the idea in [6] can be realized by the singlino-like neutralino dark matter (hereafter the 


singlino-like neutralino is simply called singlino): 


(i) The singlino dark matter annihilates to the light singlet Higgs bosons and the relic 
density can be naturally obtained from the interaction between singlino and singlet 


Higgs bosons; 


(ii) The singlet Higgs bosons, not related to electroweak symmetry breaking, can be light 


enough to be kinematically allowed to decay dominantly into muons or electrons 


through the tiny mixings with the Higgs doublets; 


(iii) The Sommerfeld enhancement needed in the dark matter annihilation for the expla- 


nation of Pamela result can be induced by the light singlet Higgs boson (h). 


Such an explanation of dark matter requires that the singlet Higgs field has very small 
mixing with the Higgs doublets, which implies that the singlino dark matter may remain 
hidden and irrelevant to the LHC experiments. However, we note that the singlet extension 
of the MSSM has a quite large parameter space and thus the coupling of the light singlet 
Higgs (h, a) with the doublet Higgs (the lightest one is called hsm) may be enhanced by 
other parameters. For example, through the soft term ApS? (S is the singlet Higgs field) 
with a large A,, a pair of singlet Higgs bosons may sizably couple to a doublet Higgs boson 
although the mixing between the singlet and doublet Higgs fields is small. Therefore, this 
model may allow for exotic Higgs phenomenology at the LHC. 

In this work we study the parameter space allowed by the explanation of Pamela result 
plus relic density via Sommerfeld enhancement and also consider the constraints from the 
LEP experiments. We find that although the light singlet Higgs bosons have small mixings 
with the Higgs doublets, their couplings with the SM-like Higgs boson (hsm) can be enhanced 
by the soft parameter A, and, in order to meet the stringent LEP constraints, the hsm 
tends to decay into the singlet Higgs pairs aa or hh instead of bb. This implies that the Asm 
produced at the LHC will give a multi-muon signal, hsm —> aa — 4p or hsm > hh — 4a — 
8u. 

This work is organized as follows. In Sec. II we discuss the Higgs and neutralino sectors 
in the singlet extension of the MSSM. In Sec. III we scan the parameter space allowed by 
the dark matter explanation and LEP experiments, and discuss the implication on Higgs 


phenomenology. Finally, a summery is given in Sec. IV. 


II. HIGGS AND NEUTRALINOS IN SINGLET EXTENTION OF MSSM 


The Higgs superpotential in the general singlet extension of the MSSM is given by [10] 
FP ae dte «ai te Mr dd MG 
W = pHa Ha + MSHu- Ha +08 + 7 1,5? + 3658, (1) 


where 5 is the singlet Higgs superfield while H, and H, are the doublet Higgs superfields. 
The Higgs scalar potential consists of the D-term, the F-term and the soft SUSY-breaking 
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term. Since $ is a singlet, the D-term is same as in the MSSM. The F-term from the 


superpotential is given by 
Vp = |u + AS (HU + |a|?) + [n + 1188 + AH, + Ha + 6 S?)?. (2) 
The soft SUSY-breaking terms are given by 


Voor = Min, Hu? + må, Hal? + må| S|? 


1 1 
+(BuH, : Ha + Ay Hy HaS + CNS + = BypisS” + 3 KA, S +h.c.). (3) 
So the Higgs potential reads 


V= |p AS? (|H? + Ha?) + |AHu Ha + 68°]? 


1 1 * —* 
PG 9 (|Hul? = |Hal?)? + 5 92\ Hr Ha FEH |? 


+m}, Hal? +m, |Hal? + må] S|? 


+HBuH, i Ha + AA) Ai, : AGS + àusHu å Has* 


1 1 
+CnS + z Buss" + q KA, S3 + Kus S?S* + hic.) (4) 


where g? = (g? + g5)/2 with gı and gp being respectively the coupling constant of SU(2) and 
U(1) in the SM. 
After the Higgs fields develop the vevs hu, hg and s, i.e., 


Hå uR +t ul dR tî ar Fat RT 
v2 v2 v2 


we obtain a 3 x 3 mass matrix Mp for CP-even Higgs bosons, a 3 x 3 mass matrix Ma for 


HÌ = ha + (5) 


CP-odd Higgs bosons and a 2 x 2 mass matrix Me for the charged Higgs bosons: 


(1) The CP-even Higgs mass matrix in the basis (Hur, Har, Sr) is given by 


Mari = gh + cot 8 [A8(Ay + Ks + us) + Bul, (6) 
Maz = g*hå+ tan 8 [As(Ax + Ks + us) + Bul, (7) 
Mn,33 = A(Ay + Us) ul — aE (hi, + h2) ks(A + 4ks +35) = - (8) 
Mn 12 a (2d? = 9’ )huha = As(A) ag + Ls) — Bu, (9) 
Mraz = 2A(u + As)hy — Ma(Ayx + 2K8 + ps), (10) 
Mno3 = 2A(u + As)ha — Ahu (Ax + 2K8 + ps), (11) 
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where tan 9 = hu/ha. This mass matrix can be diagonalized by a rotation 


hy Hur 
ha | =U | Har (12) 
hg SR 


with an orthogonal matrix U. The mass eigenstates are ordered as Mp, < Mhs < Mha. 
In the MSSM limit (À, 7, Hs, k 2 0 and h3 ~ Sp) the elements of the first 2 x 2 
sub-matrix of U are related to the MSSM angle a as 


U= cosa, Uz, = sina , 


Up = -sina , Uo, = cosa . (13) 


The CP-odd Higgs mass matrix Ma in the basis (Hur, Har, Sr) is given by 


Maj = cot B[AS(Ax + Ks + us) + Bul, (14) 
Ma2 = tan B[AS(Ay + Ks + us) + Bul, (15) 
huha 


Mags = 4AKhyha + A(Ay + Us) 


S 


-AEA + DÅ) — roA + ta) = Å — 2Bate (16) 
Ma,2 = A8(A, + KS + us) + Bu, (17) 
Mais = Aha(Ay — 268 — ps), (18) 
ud) (19) 


The diagonalization of this mass matrix can be performed in two steps. The first step 


is to rotates into a basis (A, G, Sr) with G being a massless Goldstone mode: 


Aut cos 3 —sin 8 0 Å 
Har | = | sing cosB 0 CAR (20) 
Sr 0 0 1 Si 


Dropping the Goldstone mode, the remaining 2 x 2 mass matrix in the basis (Å, Sr) 


is given by 
Mati = (tan 8 + cot 3)[As(A, + Ks + us) + Bul, (21) 
Man = 4Arhuha + MAx + po) 202 — AEC + 3) 
—5(3Aq + pu) = Cl — 2Bate (22) 


Ma2 = M/ h? + h? (Ax — 2K8 — ls). (23) 
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SS 


SS 


It can be diagonalized by an orthogonal 2 x 2 matrix P’ and the physical CP-odd 


states a; are given by (ordered as ma, < Ma,) 


ay = PÅ Dg Pad = PÅ, (cos BHur + sin GHar) a Po (24) 


a2 = PĀ + PS1 = P, (cos BHaur + sin Har) + Pr (25) 


The charged Higgs mass matrix M+ in the basis (H as Hi) is given by 


2 t 1 
M: = (astas HKS + flg) + But huh — 5) cot 8 , (26) 


which gives one eigenstate H* of mass TrM4+ and one massless goldstone mode G*: 


H= 


u 


cos BH=* — sin GF , 


Hj = sin 8H* + cos BG . (27) 


The neutralino mass matrix Mo can be read from the Lagrangian 
1 1 343 
Lis 5 MAA + 5 MMM 
1 
tuputa + Asbuda + hupat + haput) — (K8 + 5 Hs) Vss 


er XM (hud — had) — i A (huh? — had), (28) 


where A; is the U(1)y gaugino and \3 is the neutral SU(2) gaugino. In the basis 
y? = (id, ida, VI, YI, Ys) we obtains 


L= -5 V' Mol(Y')T + h.c., (29) 
where 
M 0 2 a 0 
M,-22 > 278 0 
Mo = 0 Ne) —àha |- (30) 
0 Nie 
2KS + [Lg 


Diagonalizing this mass matrix, one obtains 5 mass eigenstates (ordered in mass) 


XP = Nyt}. (31) 


III. EXPLANATION OF PAMELA AND IMPLICATION ON HIGGS DECAYS 


In our study the lightest CP-odd neutral Higgs boson a is singlet-dominant, while for the 
CP-even neutral Higgs bosons the lightest one hj is singlet-dominant and the next-to-lightest 


hə is doublet-dominant. We use the notation: 
a=a, h=h, hsm =ho. (32) 


As discussed in [10], when the lightest neutralino X? in Eq.(31) is singlino-dominant, it can 
be a perfect candidate for the dark matter. As shown in Fig.1, such singlino dark matter 
annihilates to a pair of light singlet Higgs bosons followed by the decay h — aa (h has very 
small mixing with the Higgs doublets and thus has very small couplings to the fermions). In 
order to decay dominantly into muons, a must be light enough. Further, in order to induce 
the Sommerfeld enhancement, h must also be light enough. From the superpotential term 
«5° we know that the couplings hy? X? and ay? Y? are proportional to «s. To obtain the relic 


density of the dark matter, x should be O(1). 


x? TTT 
titrs 


rhage 
FELT * 
Alitite 
ILL Le 
ILL Le 
EPES 


xf 


FIG. 1: Feynman diagrams for singlino dark matter annihilation where Sommerfeld enhancement 


is induced by exchanging h. 


Since hk, a must be singlet-dominant and ¥? must be singlino-dominant, this implies small 
mixing between singlet and doublet Higgs fields. From the superpotential in Eq.(1) we see 
that this means the mixing parameter Å must be small enough. On the other hand, the 


smallness of À is also required by the lightness of hj and a, whose masses are approximately 


given by 
huh C 
Ma 33 = KS Aas + Hs) “ — EN + hå) + (4. + 4KS + 3g — Eu , (33) 
KS KS KS 
huha H 2 2 Cy | 2Bspts 
Ma,22 ~ KS Aas ST: Hs) KE = Azza (Pu + hi) = 3A, + Hs + ee + ae i (34) 


In the following we scan over the parameter space. We modify the package NMSSM Tools 


[11] and use it in our calculations. As discussed above, Å must be small enough in order to 
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get a singlino-dominant Y? and singlet-dominant h, a (we checked from our scan that Å must 
be smaller than 0.01 in order to get Mma < 0.5 GeV and mn < 20 GeV). So in our following 
scan we fix Å = 103. Further, « is taken as 0.5, and for the squark sector the soft masses 


and the trilinear terms are fixed as 500 GeV. Other parameters vary in the ranges: 


—500 GeV < C, u, us, B, Ay, Mi, Ma < 500 GeV 
—(500 GeV)? < 7 < (500 GeV)?, s< 500 GeV, 2< tan < 40. (35) 
In order to get small M),33 and Mq22, the third terms in Eqs.(33,34), which are not sup- 


pressed by a small A, must also be small. Therefore, in our scan we require parameters Å, 


and B, to be in the ranges: 


C 
An € (års gt <1) + 20GeV, (36) 
Cn 2 
2Byptg € | —3A,KS — HKS — =e + (3GeV) (37) 


In addition, we consider the following constraints: 


(i) The constraints from the LEP experiments, which include the LEP1 bound on invisible 
Z decay and the LEP2 direct searches for Higgs bosons; 


(ii) ma, < 0.5 GeV; 


(iii) The singlino-like ¥? to give the dark matter relic density Qgoh? in the range 0.01-0.2, 


which can be calculated from the approximate formula [10] 


Py EE Å ax i (38) 
xi K 200GeV/ ` 


To calculate the Sommerfeld enhancement we follow [6] to numerically solve the Schrödinger 


equation 
1 @& k? 
"om ge +V(r)x= JMX (39) 


with the boundary condition (r — oo) 
x(r) — sin(kr + ô), (40) 


where M and k are respectively the mass and momentum of the dark matter particle. V(r) 


is the Yukawa potential induced by exchanging h and is given by 


V(r) = -— (41) 


The Sommerfeld enhancement is then given by 


2 


a (22) 


r- |20 


The survived points are displayed in different planes in Figs.2-6. We see from Fig.2 that 
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FIG. 2: The scatter plots showing the decay branching ratios a > tu (muon), a — gg (gluon) 


and a — s5 (s-quark) versus ma for Å = 103. 


in the range 2m, < Ma < 2m,, a decays dominantly into muons. From Fig.3 it is clear that 
h can be as light as a few GeV, which is light enough to induce the necessary Sommerfeld 
enhancement as shown in Fig.4. In the calculation of the Sommerfeld enhancement, we 
assumed the dark matter move with a velocity 150 km/s. 

The fit to Pamela result has been given in [10]. As shown in Table I in [10], for the 
parameter space in Figs.2-4 with 2m, < Ma < 2m, and mn as light as a few GeV (so 
the Sommerfeld enhancement factor is large enough), the Pamela positron excess can be 
naturally explained. 

In Fig.5 we show the branching ratios of hsm decays. We see that in the allowed parameter 
space ham tends to decay into aa or hh instead of bb. This can be understood as following. 
The MSSM parameter space is stringently constrained by the LEP experiments if hsm is 
relatively light and decays dominantly to bb, and to escape such stringent constraints hsm 


tends to have exotic decays into aa or hh. As a result, the allowed parameter space tends to 
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FIG. 3: Same as Fig.2, but showing Mp and mMpg,/ Versus Ma. 
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FIG. 4: Same as Fig.2, but showing the Sommerfeld enhancement factor induced by h. 


favor a large A, as shown in Fig.6, which greatly enhances the couplings Asmaa and hgyhh 
through the soft term «A,,.S° although S has a small mixing with the doublet Higgs bosons. 
Such an enhancement can be easily seen. Take the coupling Asmhh as an example. The soft 
term KA,,S° gives a term kA, 57% which then gives the interaction KA, U? U23 hgwhh because 


SR = Ush + Uz3ha + Uz3h3 with hy =hand ho = hsm (see Eqs.12 and 32). Although the 
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FIG. 5: Same as Fig.2, but showing the branching ratio of hsm decays. The ’o’ (blue), °x’ (green) 


and ’+’ (red) denote the branching ratios of hsm — aa, hsm — hh and hsm > bb, respectively. 
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Br(Agyaa, hh) 
FIG. 6: Same as Fig.2, but showing |A,,| versus the branching ratio of hsm > aa, hh. 


mixing U?,U 3 is small for a small A, a large A, can enhance the coupling hsmhh. 
The SM-like Higgs boson hgm will be intensively searched at the LHC and its dominant 
decay mode in the MSSM is bb. In the singlet extension of the MSSM, its dominant decay 


mode may be changed to aa or hh, as shown in our above results. Such new decay modes 
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will give a multi-muon signal for hsm at the LHC, i.e., hsm —> aa — 4u or hsm > hh — 


4a — 8u. So the phenomenology of hsm will be quite different from the MSSM predictions. 


Finally, we make some remarks regarding our results: 


(1) 


The recent DO search for h — aa — 4u or 2427 channel obtained null results, which 
constrained the parameter space for the CP-odd Higgs a in the mass range of 3.6-9.5 
GeV [12]. But they do not constrain the parameter space considered in our analysis 
because we considered a much lighter CP-odd Higgs a with a mass below 0.5 GeV. 
Also, as pointed in [10], such a light a is allowed by T(3s) — ya —> yutu” [13] and 


K* — nta — at utp [14] because in our scenario a is over dominated by singlet. 


In the allowed parameter space displayed in our results, the mass of the SM-like Higgs 
boson hsm is rather below its theoretical upper bound (about 135 GeV in the MSSM). 
The reason is that, in order to push up its mass, the loop effects of heavy stops are 
needed ( note that in the singlet extension the tree-level upper bound can be enhanced 
by a term proportional to A, which is very small in our scenario). In our calculations 
the soft mass parameters in the squark sector are fixed to be 500 GeV and hence the 
stops are not heavy enough to push the mass of hsm up to 135 GeV. Of course, we can 
choose heavy stops to push up the mass of Asm, in which case the allowed parameter 
space displayed in our results (with a relatively light hsm decaying dominantly into 


aa or hh) can still survive. 


For the specified singlet extensions like nMSSM and NMSSM [15], the explanation 
of Pamela and relic density through Sommerfeld enhancement is not possible. The 
reason is that the parameter space of such models is stringently constrained by various 
experiments and dark matter relic density [16], and, as a result, the neutralino dark 
matter may explain either the relic density or Pamela, but impossible to explain both 
via Sommerfeld enhancement [17]. For example, in the nMSSM various experiments 
and dark matter relic density constrain the neutralino dark matter particle in a narrow 


mass range [16], which is too light to explain Pamela. 
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IV. SUMMARY 


The singlino-like dark matter in the MSSM extended by a singlet Higgs superfield can 
give a perfect explanation for both the relic density and the Pamela result through the 
Sommerfeld-enhanced annihilation into singlet Higgs bosons (a or h followed by h — aa) with 
a being light enough to decay dominantly to muons. In this work we analyzed the parameter 
space allowed by such a dark matter explanation and also considered the constraints from 
the LEP experiments. We found that although the light singlet Higgs bosons have small 
mixings with the Higgs doublets in the allowed parameter space, their couplings with the 
SM-like Higgs boson hsm can be enhanced by the soft parameter A, and, in order to 
meet the stringent LEP constraints, the hsm tends to decay into the singlet Higgs pairs 
aa or hh instead of bb, which will give a multi-muon signal for hsm produced at the LHC, 


hsm > aa — 4u or hsm > hh — 4a — Bu. 
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